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' ABSTRACT 
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Aims. The aim of this work is to identify the so far unknown low mass stellar population of the ~ 2Myr old Cygnus OB2 star forming region, 
and to investigate the X-ray and near-IR stellar properties of its members. 

Methods. We analyzed a 97.7 ksec Chandra ACIS-I observation pointed at the core of the Cygnus OB2 region. Sources were detected using 
the PWDetect code and were then positionally correlated with optical and near-IR catalogs from the literature. Source events were then ex- 
] traded with the Acis Extract package. X-ray variability was characterized through the Kolmogorov-Smirnov test and spectra were fitted using 
• absorbed thermal plasma models. 

r- H Results. We detected 1003 X-ray sources. Of these, 775 have near-IR counterparts and are expected to be almost all associated with 
\^ , Cygnus OB2 members. From near-IR color-color and color-magnitude diagrams we estimate a typical absorption toward Cygnus OB2 of 
■ A v « 7.0 mag. Although the region is young, very few stars (~4.4%) show disk-induced excesses in the near-IR. X-ray variability is detected in 
~13% of the sources, but this fraction increases, up to 50%, with increasing source statistics. Flares account for at least 60% of the variability. 
. Despite being generally bright, all but 2 of the 26 detected O- that early B-type stars are not significantly variable. Typical X-ray spectral 
I ■ parameters are log N H ~ 22.25 (cnr 2 ) and kT ~ 1.35 keV with lcr dispersion of 0.2 dex and 0.4 keV, respectively. Variable and flaring sources 
O | have harder spectra with median kT=3.3 and 3.8 keV, respectively. OB stars are typically softer (kT~0.75 keV). X-ray luminosities range 
+3 ■ between 10 30 and 10 31 erg s -1 for intermediate- and low-mass stars, and 2.5xl0 30 and between 6.3xl0 33 erg s -1 for OB stars. 
^ ' Conclusions. The Cygnus OB2 region has a very rich population of low-mass X-ray emitting stars. Circumstellar disks seem to be very scarce. 
. . , X-ray variability is essentially related to the magnetic activity of low-mass stars (M/M o ~0.5 to 3.0) which display X-ray activity levels compa- 
^ ■ rable to those of Orion Nebular Cluster (ONC) sources in the same mass range. 
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1. Introduction 

The Cygnus OB 2 association is one of the richest star form- 
ing regions in the Galaxy, containing a large population of O- 
and B-type stars, some of which are among the most mas- 
sive stars known jTorres-Dodeen et"ail Il99ll) . For its rich- 
ness Cygnus OB2 has been in the pa st considered a "young 
globular cluster" (Knodlseder, 2000), although unc ertainties 
remai n on the real size of total cluster popu lation (IHansonl 
I2003I) . At a distance of 1.45 kpc (DM= 10.80. IHansonl 12003^ 
Cygnus OB2 lies behind the Great Cygnus Rift and is af- 
fected by large and non-uniform extinction, A v : 5 to 15. Its 
distance and absorption have so far hindered optical stud- 
ies of intermediate- and low- mass stars in the region. The 
position of known massive OB stars in the HR diagram of 
Cygnus OB2 suggests an age betwee n 1 and 3 Myr, with ~ 2 
Myr as the most probab le value jMassev & Thompson! 1 1 99 it 
iKnodlseder etaill2002h . 

Send offprint requests to: J.F.A.C. - email: facundo@astropa.unipa.it 



The Cygnus OB2 region played a remarkable role in the 
history of stellar X-ray astronomy as an early X-ray obser- 
vatio n, performed with the Einstein satellite and analyzed 
by iHarnden et al first revealed that early-type stars 

are intense X-ray emitters. This work was also the first to 
suggest the idea that X-rays in OB stars are produced in 
shocks within stell ar winds, the start ing point for the X-ray 
emission model of iLucv & W hite ( 1980). Subsequent ROSAT, 
ASCA, and Chandra-HETG observations of massive stars in 
the Cygnus OB2 region have improved our understanding of 
the mechanis ms responsible for t he X-ray emission of O- type 
stars JWaldron et alll2004fl99llKitamoto & Mukai[ |l99fj) . 

These latter studies of the region have focused on the 
four most known OB-type stars of the region, Cyg#5, 8, 9, 
and 12, but have left the X-ray properties of the other high- 
, intermediate-, and low-mass stars in the region unexplored. 
Indeed the intermediate/low-mass population has not even been 
identified yet. Because of the much higher X-ray luminosity of 
young stars (e.g. log(L x ) ~ 30- 31 erg/s, for M = 1.0- 3.0 M Q ) 
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with respect to older field stars (e.g. lGetman et alll2005al) . X- 
ray observations have proved an efficient tool to select likely 
members of young Star Forming Regions (SFR). Thanks to 
its high sensitivity and spatial resolution, the Advanced CCD 
Imaging Spectrometer (ACIS) camera on board of the Chandra 
satellite ( Gar mire et all l2003_ is particularly suited for this 
kind of studies in crowded stellar fields such as Cygnus OB2 . 

In this paper we present the analysis of a deep (~ 97.7 ks) 
Chandra ACIS observation of Cygnus OB2 . In §2 we intro- 
duce our X-ray observation and describe our data calibration 
and reduction procedure. In section 3 we describe source de- 
tection, photon extraction and first characterization of X-ray 
spectra via hardness ratios. In §4 we characterize our sources 
by cross-identifying them with available optical and near-IR 
catalogs. Section 5 deals with X-ray variability while in §6 the 
results of the X-ray spectral analysis is presented. In §7 we dis- 
cuss out results concerning the X-ray characteristics of high-, 
intermediate- and low-mass stars of the region. Finally, in §8 
we summarize our results and draw our conclusions. 

2. The X-ray observation 

Cygnus OB 2 was observed with the ACIS d etector on board 
the C handra X-ray Observatory (CXO) ( Weis skopf et all 
l2002_ on 2004 January 16 (Obs.Id. 4511; PI of the observa- 
tion: E. Flaccomio). The total exposure time was 97.7 ksec. 
The data were acquired in very faint mode, to ease filtering 
of background events, with six CCD turned on, the four com- 
prising the ACIS-I array [0,1,2,3], plus CCDs 6 and 7, part 
of ACIS-S. However, these latter two CCDs will not be used 
in the following because of the much degraded point spread 
function (PSF) and reduced effective area resulting from their 
large distance from the optical axis. The ACIS-I 17'xl7' field 
of view (FOV) is covered by 4 chips each with 1024x1024 
pixels (scale 0.49" px _1 ). The observation was pointed toward 
R.A.=20 h 33 m 12.2 s and DEC=+41° 15' 00.7", chosen so to 
maximize the number of stars in the FOV and to keep two 
source-dense regions close to the optical axis, where the PSF is 
sharper. Figure^shows Cygnus OB2 as seen in X-rays by our 
ACIS-I observation. 

2.1. Data reduction 

Data reduction, starting with the Level 1 event list provided 
by the pipeline processing at the CXO, was performed using 
CIAO 3.2.2 1 and the CALDB 3.1.0 set of calibration files. We 
produced a level 2 event file using the acis_process_event CIAO 
task, taking advantage of the VF-mode enhanced background 
filtering, and retaining only events with grades = 0,2,3,4,6 and 
status =0. Photon energies were corrected for the time depen- 
dence of the energy gain using the corr_tgain CIAO task. 
Intervals of background flaring were searched for, but none 
were found. We will hereafter assume a constant background. 
To improve sensitivity to faint sources, given the spectrum of 
the background and that of typical sources, we filtered out 
events outside the [500:8000] eV energy band. 

1 http://cxc.harvard.edu/ciao/ 



An exposure map, needed by the source detection algorithm 
and to renormalizes source count-rates, was calculated with 
the CIAO tool mkexpmap assuming a monochromatic spectrum 
(kT = 2.0keV), 2 . 

3. Analysis 

In this section we describe the first steps taken for the analy- 
sis of the ACIS-I data. We discuss the detection of X-ray point 
sources ( 33. 1> . the definition of source and background event 
lists ( 33. 2> . and a first spectral characterization of source spec- 
tra through hardness ratios ( 33.3> . 

3.1. X-ray Source Detection 

Source detection was perfor med with the Palermo Wavelet 
Detection code, PWDetect 3 (Damia nTetail Il997bl) . on the 
level 2 event list restricted to the [500:8000] eV energy band. 
PWDetect analyzes the data at different spatial scales, allow- 
ing the detection of both point-like and moderately extended 
sources, and efficiently resolving close sources pairs. The most 
important input parameter is the detection threshold (SNR), 
which we establish from the relationship between threshold, 
background level of the observation, and expected number of 
spurious detections due to Poisson noise, as determined from 
extens ive simulations of source-free fields (cf. Damiani et all 
Il997al) . The background level was determined with the back- 
ground command in the ximage 4 package. The method amounts 
to dividing the image into equal-size boxes, discarding those 
that, according to several statistical criteria, are contaminated 
by sources, and finally computing the mean level of the remain- 
ing ones. We obtain that our ACIS-I observation comprises a 
total of about 1.58xl0 5 background photons. This background 
level translates into a SNR threshold of 5.2 if we decide to ac- 
cept one spurious detection in the FOV, or into SNR > 4.5 if we 
decide to accept 10 spurious detections. The first choice results 
in the detection of 868 sources and the second one in 1054 
sources. By accepting an extra ~9 spurious detections in the 
FOV we thus gain about 177 new reliable sources. Considering 
moreover that 10 spurious sources amounts to only ~1% of the 
total number of detections, we decided to adopt the second less 
conservative criterion. 

After a careful visual inspection of the initial source list, we 
rejected a total of 51 detections which we considered spurious: 
39 were produced by different instrumental artifacts (i.e. out- 
of-time events, CCD gaps, detector edges, etc) that were not in- 
cluded in the simulations of source-free fields used to establish 
the SNR threshold, but that are easily recognized. The remain- 
ing 12 were multiple detections of the same sources with dif- 
ferent spatial scales. 5 In total, our final list of X-ray sources in 

2 http://asc.harvard.edu/ciao/download/doc/expmap_intro.ps 

3 See http://www.astropa.unipa.it/progetti_ricerca/PWDetect 

4 http://heasarc.nasa.gov/xanadu/ximage/ximage.html 

5 This effect occurred for only 3 sources at large off-axis angles, 
where the PSF is particularly elongated and thus significantly differ- 
ent from the symmetric PSF assumed by PWDetect. Note that, quite 
obviously, this effect is also not included in the simulations of source- 
free fields. 
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Fig. 1. Upper panel: color-coded ACIS-I image of the 17' x 17' field in Cygnus OB2 (see color version in the electronic edition). 
Kernel smoothing was applied to highlight point sources. Energy bands for the RGB image are [0.5:1.7], [1.7:2.8], and [2.8:8.0] 
ke V for the red, green, and blue colors, respectively. The four panels at the bottom show four 30" x 20" areas of the ACIS field at 
four different distances from the optical axis: from left to right: 0.7, 4.7, 5.4, and 7 arcmin. They illustrate the dependence of the 
PSF quality on off-axis and show (black continuous lines) the source photon extraction regions established with ACIS Extract. 
Figures next to these contours indicate the source identification numbers (from Table 1) and the encircled energy fractions. 



the Cygnus OB2 region contains 1003 X-ray detections, 99% 
of which are expected to be real. The first 7 columns of Tabled 
list, for each source: a running source number, name (accord- 
ing to CXC naming convention 6 ), sky position (R.A. and Dec.) 
with relative uncertainty, off-axis angle (9), significance of the 
detection (Sig.). 



6 http://cxc.harvard.edu/cdo/scipubs.html 



3.2. Photon Extraction 

For further analysis of the detected sou rces, we used ACIS 
Extract 7 (AE) v3.79 jBroos et alll2002h . an IDL based pack- 
age that makes use of TARA 8 , CIAO and FTOOLS 9 . 



7 http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 

8 http://www.astro.psu.edu/xray/docs/TARA/ 

9 http://heasarc.gsfc.nasa.gov/docs/software/ftools/ 
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The optimal extraction of point source photons is compli- 
cated by the non-Gaussian shape of the PSF and by its strong 
non-uniformity across the ACIS field of view. In particular the 
width and asymmetry of the PSF depend significantly on the 
off-axis distance (ff). The PSF is narrow and nearly circular in 
the inner fl<5' but becomes rapidly broader and more asym- 
metries at larger off-axis, as demonstrated in the bottom panel 
of Fig. [0 where we show sources detected in our data at four 
different off-axis angles. 

Extraction from circular regions containing a large frac- 
tion (e.g. 99%) of the PSF would guarantee to collect almost 
all source photons. However, given the extended wings of the 
PSF, very large regions would be needed, incurring in the risk 
of contamination from nearby sources, especially in crowded 
fields like ours. Moreover the resulting inclusion of a large 
number of background events would reduce the signal to noise 
of weak sources. On the other hand, extraction from regions 
that are too small may reduce the photon statistic for further 
spectral and timing analysis. 

In order to tackle this problem, AE begins with calculat- 
ing the shape of the model PSF at each source position using 
the CIAO task mkpsf. It then refines the initial source posi- 
tions (in our case estimated by PWDetect assuming a symmet- 
ric PSF) by correlating the source images with the model PSFs. 
Following AE science hints 10 , this last procedure was only used 
for those sources lying at off-axis larger than 5 arcmin (464 
sources), while for the rest of source (539 sources) we sim- 
ply adopt data-mean positions. Coordinates listed in Tabled 
as well as their lcr uncertainties, are the result of this pro- 
cess. We verified that these positions are an improvement over 
those computed by PWDetect, especially at large off-axis an- 
gles, comparing the offsets between X-ray sources and coun- 
terparts in the 2MASS catalog (see §4.2). 

After re-computing positions, AE defines source extrac- 
tion regions as polygonal contours of the model PSF contain- 
ing a specified fraction of source events (fpsF)- Generally, we 
chose fpsF=90%, and computed the contours from the PSF for 
a mono-energetic source with E = 1.49 keV. For ~8% of the 
sources in the denser parts of the Cygnus OB2 field this frac- 
tion was reduced so to avoid contamination with other nearby 
sources, in the most extreme cases down to fpsF~40%. The four 
panels at the bottom of Fig.^show examples of extraction re- 
gions. Otherwise, we have detected only six sources: Id.# 17, 
60, 488, 544, 568, 729 that suffers of a pile-up fraction (f p u e -up) 
of3.96%, 17.7%, 11.6%, 2.39%, 15.6% and 4.6%, respectively. 
In these cases, we extracted events in annular regions that ex- 
clude the PSF cores for all sources that exhibits a f p ii e - U p greater 
than 2%. The inner radii was fixed at 2 arcsec 11 , while the outer 
circles were chooses as the smallest that inscribe the fpsF=99% 
contours. 

Although the ACIS-I instrumental background level is spa- 
tially quite uniform, the actual observed background varies 
substantially across the crowded Cygnus OB2 field due to the 
extended PSF wings of bright sources and to their readout 
trails. Background was therefore estimated locally for each 

10 http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide/node35 

11 http://www.astro.msfc.nasa.gov/Ch4/Ch4_15-03_Tsujimoto.pdf 



source, adopting once again the automated procedure imple- 
mented in AE, which defines background extraction regions 
as circular annuli with inner radii 1 . 1 times the maximum dis- 
tance between the source and the 99% PSF contour, and outer 
radii defined so that the regions contains more than 100 "back- 
ground" events. In order to exclude contamination of the re- 
gions by nearby sources, background events are defined from 
a "swiss cheese" image that excludes events within the inner 
annuli radii of all the 1003 sources. 

Results of the photon extraction procedure are listed in 
columns 8-16 of Tabled We give: the source extraction area 
(column 8); the PSF fraction within the extraction area, assum- 
ing E= 1 .49 keV (9); the background-corrected extracted source 
counts in the 0.5-8.0 keV band (10); the exposure time (11), the 
count rates (CR) in four spectral bands computed as the ratio 
between the source photons (corrected for fps f) and the expo- 
sure time (12-15); the median photon energy (E x ), (16). 

In summary, our 1003 X-ray detections span a wide count 
range, from 4 to ~15000 photons. Most sources are faint (e.g. 
42% have less than 20 photons). The lower envelope of the 
counts vs. off-axis plot (not shown) indicates that the minimum 
number of photons necessary for detection ranges from 4 on 
axis to ~20 close to the detector corners (0=10 arcmin). 

3.3. X-ray hardness ratios 

In order to characterize the X-ray spectra of low-statistic 
sources it is common practice to use the ratios of source counts 
in different spectral bands, i.e. X-ray hardness ratios (XHR), or 
the logar ithm of these values which can be consid ered "X-ray 
colors" .Schulz et allll989_|Prestwich et allEool . 

We use this method, dividing the full energy range into 
three bands: Soft (S x : 0.5-1.7 keV), Medium (M x : 1.7-2.8 
keV) and Hard (H x : 2.8-8.0 keV). Figure [2] shows, separately 
for sources in three detected counts ranges, the "Hard X-ray 
color" (log[M x /// x ]) vs. the "Soft X-ray color" (log[S x /M x ]). 
For reference we also plot the predicted loci for absorbed 
thermal sources with plasma temperatures between 0.5 and 
8.0 keV and N H between 10 20 and 10 23 cm" 2 . The grid 
was calculated with the Portable Interactive Multi-Mission 
Simula tor (PIMMS 12 ) for a Ray mond-Smith (RS) emission 
model dRavmond & Smifhlll977t . 

A comparison of the three panels in Fig. [2] shows that the 
position of sources with respect to the kT-Ne grid is signifi- 
cantly affected, other than by the source spectra, by the source 
statistic. The positions of individual sources with less than 20 
counts (left panel) are, for example, considerably spread-out 
because of large statistical uncertainties. Sources with higher 
statistics (20 < NetCnts < 50, central panel) are more concen- 
trated around kT ~ 2.0 keV and Nh ~ 2.3xl0 22 cur 2 . Sources 
with more than 50 photons (right panel) are even more con- 
centrated around these same values of kT and Nh. A separate 
group of soft-spectrum sources becomes however visible in the 
upper part of the grid, with typical kT~0.8 keV. Most of these 
sources are associated with known early type stars, indicated 



12 http://heasarc.gsfc.nasa.gov/docs/software/tools/pimms.html 



Table 1. Cygnus OB2 X-ray source catalog: first 35 rows. The complete table, containing 1003 rows is available in the electronic edition of A&A. 
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3 
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4 
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14 
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3.18 
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5 
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0.76 
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55 
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6 
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7 
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8 
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33 
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0.078 
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9 
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20:32:29.26 


41: 
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10.17 


6.70 
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0.82 


64 


92.6690 
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0.364 


0.202 
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2.02 


-0.44 


10 


203229.84+411454.2 


20:32:29.85 


41: 


: 14:54.27 


0.50 


7.97 


5.56 


525 


0.90 


50 


93.3914 


0.605 


0.249 


0.185 
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2.24 


-0.67 


11 


203229.91+411056.0 


20:32:29.92 


41: 


: 10:56.06 


0.72 


8.94 


2.77 


751 


0.90 


19 
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0.227 


0.080 


0.058 
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2.30 


-0.97 


12 


203230.63+410829.1 


20:32:30.64 


41: 


:08:29.11 


0.58 


10.19 


8.52 


1205 


0.90 


99 


94.5473 


1.173 
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0.478 


0.375 


2.14 


-0.08 


13 
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20:32:30.63 


41: 


:08:54.82 


0.66 


9.92 


5.21 
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0.83 


44 


96.4507 


0.552 


0.067 
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2.93 
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14 


203230.88+411029.5 


20:32:30.89 


41: 
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0.46 


8.99 


7.08 
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0.90 


75 


96.5167 


0.867 
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0.176 
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-0.07 


15 


203231.39+410955.9 


20:32:31.40 


41: 
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0.57 


9.21 


6.16 
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0.90 


60 


94.7546 


0.714 


0.259 
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2.27 


-0.01 


16 


203231.42+411335.0 


20:32:31.43 


41. 
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7.80 


3.85 
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0.90 


30 
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0.098 
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3.01 
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17 
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41 
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0.89 


3173 
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1.85 
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18 
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41: 
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0.73 
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2.01 
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0.90 


12 
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3.10 


-0.04 


19 
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0.56 


7.50 


3.55 
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0.89 


25 


94.0133 
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3.23 
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20 
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0.43 
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49 
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21 
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18 
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1.15 


-0.44 


22 


203233.47+411218.0 
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41 
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0.49 


7.77 


4.49 
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0.90 


34 


91.7675 


0.415 


0.177 


0.127 


0.111 


1.82 


-0.04 


23 


203233.53+411405.7 


20:32:33.54 


41: 


: 14:05.74 


0.38 


7.33 


5.47 


364 


0.90 


50 


95.7566 


0.593 


0.199 


0.201 


0.193 


1.92 


-0.49 


24 


203234.33+411506.3 


20:32:34.34 


41: 


:15:06.33 


0.54 


7.12 


3.71 


350 


0.90 


26 


94.3211 


0.311 


0.074 


0.056 


0.181 


3.52 


-0.44 


25 


203234.45+411154.7 


20:32:34.45 


41: 


:11:54.74 


0.58 


7.75 


3.96 


446 


0.90 


27 


94.8080 


0.326 


0.061 


0.203 


0.062 


2.06 


-0.26 


26 


203235.17+411143.4 


20:32:35.18 


41: 


: 11:43.48 


0.47 


7.70 


4.97 


430 


0.90 


39 


91.6010 


0.475 


0.196 


0.159 


0.120 


1.96 


-0.03 


27 


203235.65+41 1509.5 


20:32:35.65 


41: 


:15:09.50 


0.33 


6.88 


7.09 


312 


0.90 


70 


96.3691 


0.815 


0.228 


0.368 


0.219 


2.10 


-1.17 


28 


203235.86+411248.4 


20:32:35.87 


41: 


: 12:48.42 


0.37 


7.18 


5.22 


340 


0.90 


43 


93.8249 


0.519 


0.083 


0.289 


0.147 


2.04 


-2.14 t 


29 


203235.96+412135.6 


20:32:35.97 


41: 


:21:35.64 


0.78 


9.47 


4.51 


669 


0.83 


33 


95.8351 


0.413 


0.190 


0.064 


0.159 


1.85 


-0.11 


30 


203236.32+411900.7 


20:32:36.32 


41: 


: 19:00.76 


0.47 


7.84 


4.96 


489 


0.90 


41 


94.6415 


0.484 


0.209 


0.129 


0.146 


1.76 


-0.44 


31 


203236.38+411831.0 


20:32:36.39 


41: 


:18:31.08 


0.62 


7.59 


3.81 


403 


0.89 


26 


90.0839 


0.333 


0.242 


0.060 


0.032 


1.41 


-1.46 


32 


203236.50+411810.3 


20:32:36.50 


41. 


:18:10.36 


0.67 


7.42 


3.78 


376 


0.89 


26 


91.6356 


0.320 


0.115 


0.112 


0.093 


1.90 


-1.55 


33 


203236.62+412214.5 


20:32:36.63 


41: 


:22: 14.52 


0.26 


9.85 


14.02 


689 


0.81 


232 


97.5187 


2.949 


1.112 


0.835 


1.002 


1.93 


-4.00 t 


34 


203236.86+411944.6 


20:32:36.86 


41: 


: 19:44.62 


0.40 


8.16 


7.22 


544 


0.90 


74 


97.0381 


0.859 


0.295 


0.288 


0.277 


1.86 


-0.83 


35 


203236.84+412146.1 


20:32:36.85 


41: 


:21:46.11 


0.56 


9.48 


5.65 


660 


0.83 


46 


95.1661 


0.592 


0.278 


0.162 


0.153 


1.82 


-1.24 



Notes - t: flare-like sources. 
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Fig. 2. X-ray color-color diagrams for our ACIS sources. The three panels refer to sources in three ranges of detected counts, as 
indicated in the upper-left corners. Colors are defined for three energy bands: Soft (0.5-1.7 keV), Medium (1.7-2.8 keV) and Hard 
(2.8-8.0 keV) keV. The grids refers to predicted colors for Raymond-Smith thermal emission models with different temperatures 
(kT = 0.5, 1.0, 2.0, 4.0, and 8.0 keV) and affected by varying amounts of absorption, N H =0.01, 0.5, 1.0, 2.5, 5, and 10 xlO 22 
cm 2 . Boxes refer to known O- and early B-type stars. Note that the dispersion of points is larger for low-statistics sources. 
Typical lcr error bars are plotted in the right-bottom part of each panel. 



by squares in Fig. |2 Others, with somewhat smaller N H , are 
likely foreground stars as also argued in 84.21 

4. Optical and Near-IR counterparts 

Cygnus OB2 has been th e target of numerou s optic a l studies 
for more than 50 years J ohnson & Morg arl 1 19541 ISchultel 
Il958l) . iMassev & Thompsonl Jl99lh identified 120 candidate 
massive members in an area of about 0.35 deg 2 , on the basis 
of UBV photometry, and gave optica l spectral clas sifications 
for over 70 OB stars. More recently, Hanson (2003) published 
new spectral classification for 14 more OB candidates, improv- 
ing the massive star census of Cygnus OB2, and the estimates 
for its distance (1450pc ) and age (~2M yr). A total of 33 OB 
stars from the catalog of Hanson (2003) lie in the 0.0823 deg 2 
FOV of our X-ray observation and were thus cross identified 
with our X-ray source list (§3.2). Due to the high and variable 
absorption in the Cygnus OB2 line of sight, and to the relative 
shallowness of the available optical catalogs, the correlation 
between X-ray sources and optical catalogs is of limited use 
for the study of intermediate- and low-mass stars. We thus de- 
cided to base our characterization of the X-ray sources mostly 
on near-IR data, for which the impact of dust extinction is re- 
duced and actually comparable to that on the X-ray band. We 
made use of J (1.25 /urn), H (1.65 fim) and K s (2.17 //m) pho- 
tometry from the Two Micron All Sky Survey (2MASS) Point 
Source Catalog (PSC) 13 . 2MASS is complete to magnitudes of 
15.8, 15.1 and 14.3 mag in the J, H and K s bands, respectively. 
We limited our analysis to 2MASS sources for which the qual- 
ity flag for at least one of the three magnitudes is equal to A, B, 
C or D (cf. the 2MASS All-Sky Data Release User's Guide). 
With this restriction 1 1 sources were removed from our initial 
list of 5061 2MASS point sources in the ACIS FOV, leaving a 
total of 5050 objects. 

13 See http://www.ipac.caltech.edu/2mass 



4.1. Cross-identifications 

We began by cross identifying our X-ray source list with the 
2MASS catalog. Identification radii, Rjd, were chosen so to 
limit the number of spurious identifications due to chance 
alignments, N c h an ce, and at the same time to include a large 
number of the true physical associations, N trae . First we esti- 
mated, as a function of R^, the number of chance identifica- 
tions expected in a given sky area, A searc h, assuming (i) uncorre- 
cted NIR and X-ray positions and (ii) a uniform surface source 
density of the N2MASS 2MASS sources lying in the search area: 
N c hance(Rid)=NxAi d N 2 MASs/A S earch, where A id =7rR 2 d is the area 
of identification circles. Ntr U e(Rid) can instead be estimated 
from the observed total number of identifications at any given 
Rid as Nid(Rid)-N c hance(Rid)- We chose the identification radius 
as the largest for which N true >N c h an ce- 

Because the Chandra PSF, and therefore the position un- 
certainty of X-ray sources, depends mainly on the off-axis, we 
perform this analysis in four different ranges of off-axis: [0-2), 
[2-4), [4-7) and > 7 arcmin. Our identification radii in these 
four regions are 1.0, 1.5,2.1, and 2.7 arcsec, respectively. 

Before performing the final identifications we searched for 
possible systematic differences between the X-ray and 2MASS 
positions. We first performed a preliminary cross-identification 
and compared the coordinates of identified pairs. A small sys- 
tematic offset between the two catalogs was found. We thus 
shifted the Chandra coordinates and performed a new cross- 
identification, repeating this process iteratively until the off- 
set was reduced to ~0.01 arcsec, i.e. much smaller than the 
statistical errors on source positions. In the end the offset be- 
tween the two catalogs was: A(RA x _2mass) — +0.12"+0.10" 
and A(DEC x _2mass) = +0.49"±0.17". The result of the final 
identification is shown in Table [2 where we lists, in the first 
7 columns, the X-ray and near-IR identifiers of cross-identified 
sources, the offset between the two positions, and the J, H, K, 
magnitudes from 2MASS. A total of 775 X-ray sources out of 
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the 1003 in our list were identified with 2MASS objects. Three 
X-ray sources (#12, #148, and #449) were identified with two 
2MASS objects but, after checking the positions visually, we 
adopted the closest of the counterparts. The fraction of iden- 
tified X-ray sources appears to increase as we move to larger 
off-axis. In the four annular regions defined above these frac- 
tions are: 76/116 (65%), 227/304 (74%), 294/374 (79%) and 
178/209 (85%), in order of increasing off-axis. This trend can 
be attributed to the already mentioned J3.2t dependence of the 
ACIS-I sensitivity with off-axis angle and to the fact that more 
intense X-ray sources are more likely to have a near-IR coun- 
terpart than fainter ones. 

With respect to the expected number of chance identifica- 
tions we estimate, for off-axis ranges [0:2), [2:4), [4:7), and 
>7 arcmin, and with the formula given above, no more than 
1.6, 7.9, 16.4, and 12.5, respectively for a total of < 39. Note 
however that the assumption that X-ray and 2MASS source 
lists are fully uncorrected is quite certainly not true (other 
than for the 10 expected spurious detections), so that this num- 
ber can actually on ly be considered as a loose upper limit (cf. 
iDamiani et a"fll2003l) . 

We estimate the expected number of extragala ctic so urces 
in our detection list following iFlaccomio et alJ < 2006:|). We 
consider the ACIS count-rates of no n-stellar sources in the 
Chandra Deep Fie ld North (CDFN, Alexander et all 120031 
Barg er et all 120031) and estimate absorption corrected count- 
rates assuming Nh = 1.54xl0 22 (from A v = 7.0, see 34. 2t us- 
ing PIMMS and assuming power-law spectra with index F. We 
then compare these count rates with upper limits taken at ran- 
dom position in the ACIS FOV. For F between 1 and 2 we ob- 
tain 61 to 87 expected extragalactic sources. Given the intrin- 
sic near-IR fluxes of these sources and the absorption toward 
Cygnus OB2, they ar e expected to be amon g the 228 without 
NIR counterparts (c.f. lFlaccomio et alll2006h . 

Finally we identifie d our X-ray source list with the 33 early 
type stars listed bv lHansonl J2003I) and lying in our ACIS FOV. 
Of these, all the 20 O-type stars result associated with an X-ray 
source, while of the 13 B-type stars only 6 are detected. 

4.2. Near-IR properties of the X-ray sources 

We now investigate the NIR properties of the X-ray sources 
with a 2MASS counterparts. For this purpose we restrict our 
analysis to sources with high quality photometry, i.e. those for 
which the quality flag (see 2MASS documentation) is 'AAA', 
or for which uncertainties on the J, H and K s magnitudes are all 
lower than 0.1 mag. With these requirements the total number 
of IR sources in the ACIS FOV is reduced from 5050 to 2187. 
Counterparts of X-ray sources were selected only on the basis 
of their magnitude errors (<0.1 mag), yielding 519 sources out 
of the original 775. 

Figure [5] shows the J-H vs. H-K s color-color (CC) dia- 
gram for all the selected 2MASS objects in the ACIS FOV, 
both X-ray detected and undetected. We also plot for compari- 
son the MS dKenvon & Hartma nn[ll995|). the Classical T-Tauri 
Stars (CTTS1 locus of iMever et all dl997l) . and three redden- 
ing vectors starting from these loci and with slope (A Ks /E(H - 



K,;)=0 .125) taken from the extinction law given bv IHansoi] 
J2003I) . Cygnus OB2 members with purely photospheric emis- 
sion should lie in this reddening band. However, as the line of 
sight toward Cygnus OB2 is not far from that to the Galactic 
center, many field interlopers are also expected in the same 
band. Young stellar objects (YSOs), such as Classical TTauri 
and Herbig Ae/Be stars, because of the NIR excess emission 
originated in the inner parts of their circumstellar disks, are of- 
ten found to the right of this band, i.e. in the CTTS locus. 




0.0 0.5 1.0 1.5 

H-K s 



Fig. 3. JHK, color-color diagram of objects in the 
17'xl7'CygnusOB2 FOV with high quality 2MASS photom- 
etry. Filled circles and small points refer to X-ray detected 
and undetected 2MASS sources, respectively. Detected known 
O- and B-type stars are indicated by s quares. We also show 
for reference the main-sequence from Kenvon & Hartmann 
dl995h . the CTTS locus of IMever et all I 19971) and reddening 
vectors (dashed lines) with length corresponding to A v = 15 
mag. Note: i) the low number of X-ray sources lying in the 
CTTS loci; ii) the peculiar position of Cyg#12 (labeled), one 
of most intrinsically luminous stars in the Galaxy, indicative of 
a high reddening. 

We note here, however, that very few X-ray sources, i.e. 
likely Cygnus OB2 members, have colors consistent with the 
(reddened) CTTS locus. If we neglect the uncertainties on data 
points, 23 sources lie in the reddening band of the CTTS lo- 
cus. We do not include here the peculiar supergiants B5 Ie 
Cyg#12 and discuss its nature in 37.11 A total of 23 stars in the 
CTTS reddening band means a fraction of 23/5 19~ 4.4% with 
respect to all the X-ray sources in the CC diagram. We com- 
pare this fraction with the one observed in the ONC, adopting 
Chandra Orion Ultradeep Project (COUP) sources in the same 
mass range 14 that we reach in Cygnus OB2 (i.e. M> 1 M , cf. 

14 We only consider COUP sources for which mass estimates are 
given bv lGetman et alJ 1 2005a) 
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Table 2. Near-IR counterparts of Cygnus OB2 X-ray sources and near-IR stellar parameters: first rows. The complete version is 
available in the electronic edition of A&A. 



N x 


X-ray - 2MASS counterpart 


Off. 




2MASS photometry 




. 'I" 

At 


# 


CXOACJ+ 


2MASS J+ 


(") 


J 


H 


K s 


Ph. qual' 


(mag) 


1 


203225.09+411019.6 





— 














— 


2 


203225.55+410847.3 


20322545+4108473 


1.56 


14.76+0.06 


14.23±0.07 


13.95±0.07 


AAA 


8.54 


3 


203225.97+411054.7 


20322597+4110547 


0.24 


15.24±0.05 


13.65±0.03 


12.69±0.03 


AAA 


7.71 


4 


203227.51+411358.3 


20322756+4114001 


1.61 


16.00±0.06 


14.46±0.04 


13.90±0.06 


AAA 


0.16 


5 


203227.62+410831.4 


20322754+4108336 


2.41 


15.91±0.06 


14.39±0.03 


13.82±0.05 


AAA 


N/A 


6 


203227.70+411317.0 


20322768+4113169 


0.58 


14.96±0.03 


13.49±0.02 


12.88±0.03 


AAA 


2.21 


7 


203228.88+410807.5 
















8 


203229.12+411400.7 


20322913+4114012 


0.26 


16.74+0.16 


15.36±0.08 


14.72+0.11 


CAB 


9.01 


9 


203229.26+410850.1 


20322928+4108494 


1.03 


12.26±0.01 


11.49±0.01 


11.07±0.01 


AAA 


8.52 


10 


203229.84+411454.2 


20322985+4114536 


0.83 


17.44 


15.43+0.10 


14.37±0.09 


UAA 


N/A 


11 


203229.91+411056.0 


20322994+4110575 


1.29 


15.26±0.05 


14.04±0.04 


13.53±0.04 


AAA 


N/A 


12 


203230.63+410829.1 


20323051+4108284 


2.19 


15.47±0.07 


13.71 


13.09 


AUU 


4.19 


13 


203230.63+410854.8 


20323064+4108565 


1.49 


15.31±0.05 


13.68+0.03 


13.01±0.03 


AAA 


5.43 


14 


203230.88+411029.5 
















15 


203231.39+410955.9 


20323143+4109558 


0.58 


14.52±0.05 


13.04 


12.46 


AUU 


9.69 



t2MASS photometric quality flags for the J, H and, K, bands: "A" to "D" indicate decreasing quality of the measurements, "U" that the value 
is an upper limit. See 2MASS documentation for details. £N/A indicates stars with unconstrained A v (see text). 



50: 19/92~20.6%. The difference between the two fractions 
of stars with near-IR detected disks is significant, a factor of 
~4.7. It could be due to two different reasons: 

- Cygnus OB2 is older than the ONC (~ 2 vs. -lMyrs) 
and it s disks might have dissipated. According to Hillenbrand 
J2005h . the fraction of stars with disks detectable in the near-IR 
should decrease between 1 and 2 Myr by a factor of 2-3. This, 
together with the statistical uncertainties in the disk fractions 
might explain the observed difference. 15 

- The effect of disk photo-evaporation by the intense UV 
radiation field of the hot m assive OB stars. Th is mechanism, 
predicted theoretically (e.g. Ada ms et all 120041) an d observed, 
e.g., in the ONC proplyds (e.g. Ball v et all l2000l) . is thought 
to be effective in massive star forming regions, such as Cygnus 
OB2, having a large number of UV emitting O- and B-type 
stars. A shortening of the disk lifetime with decreasing distance 
from O- stars has been re cently observed in the massive star 
forming region NGC6611 JGuarcelloet all [2006. in prep.). 

Figure0]shows the K s vs. H-K s color magnitude (CM) dia- 
gram for the same stars plotted in Fig. [3] We also show for ref- 
erence the expected cluster locus: the intrinsic K s magnitudes 
and H-K s colors for stars earlier t han B 5V were taken from 
the M S calibration of [Knodlseder (2000) and lBessell&Bretll 
(1989), respectively. For later spectral types (masses between 
0.1 and 7 M Q ), we adopted the 2Myr isochrone from lSiess et all 
(2000), converted to the observational plane using the calibra- 
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14 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

H-K s 

Fig. 4. CM diagrams of the Cygnus OB 2 region. Symbols as 
in figure[3] The two parallel curves indicate the expected clus- 
ter loci for the assumed distance and for a mean reddening of 
Ay=0.0 and 7.0 mag. They represent the MS for (un-reddened) 
K s q < 10.3 and the 2 Myr isochrone for fainter magnitudes. 




A direct comparison with the results o flHillenbranj<2005h is how- 
ever not possible because her estimate of disk lifetime (i) refers to 
star in the 0.3-1.0 M Q , i.e. less massive than the ones we observe in 
Cygnus OB2, (ii) is based on a different, and more efficient, indica- 
tor of disks presence than the one we can use here, i.e. the H-K color 
excess measured with respect to the photospheric value as determined 
from spectral types. 



tion given by iKenvon & Hartmannl & 19951) . The adopted MS 
and 2Myr isochrone overlap satisfactorily. 

In order to estimate the typical visual absorption of cluster 
members we calculated the distance of each X-ray source to the 
cluster locus along the reddening direction. Resulting A v values 
for individual sources are listed in the last column of table|2] 
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Note that for 10.25<K s o < 1 1 .35 the absorption cannot be con- 
strained because the reddening vector intersect the cluster locus 
more than once. We computed median absorption values in two 
luminosity ranges: for K s o > 1 1 .35 we obtained A v =7.0 mag, 
while for K s o < 10.25, considering only the known OB stars, 
we obtained A v =5.63 mag. We note that the above estimations 
depend on the reliability of the assumed cluster locus and on 
the assumption that the H and K magnitudes are not signifi- 
cantly affected by disk-induced excesses. As observed above, 
given the paucity of stars with excesses, the latter appears to 
be a good approximation for Cygnus OB2 stars. It is anyway 
comforting that our estimate for the OB type stars is in fairly 
good agreement w ith the median of the extinctions computed 
bv iHansonl (12003) using photometric and spectroscopic data: 
A v =5.7 mag. The lower extinction derived for high mass stars 
with respect to lower mass ones may indicate that the strong 
winds and radiation field of massive stars may have cleared 
their surrounding environment. 

We note that in both the CC and CM diagrams, some X- 
ray sources, of the order of 15-20, lie close to the un-reddened 
cluster loci. These are likely to be foreground MS stars and 
thus to contaminate the sample of X-ray detected cluster mem- 
bers. This conclusion is corroborated by the X-ray hardness- 
ratio analysis: eight of these stars have indeed soft and rela- 
tively unabsorbed spectra (the stars to the right of the OB stars, 
in the upper-right corner of Fig. fright), as expected from fore- 
ground field stars. 

In the following ({JTJl we will correlate the X-ray properties 
of our sources with stellar parameters derived from the avail- 
able optical and near-IR data. We obtain an estimate of stellar 
masses for 682 2MASS counterparts from the J-band magni- 
tudes (limited to quality flags "A" to "D") and the mass vs. J 
mag. relationship appropriate to the cluster age, distance and 
extinction, obtained as described above for the cluster locus in 
the H, H-K diagram. 16 The relation is degenerate in two ranges 
of J, corresponding to O.2-O.4M and 2.8-5.3M Q . Four and 28 
X-ray detected stars lie in these ranges, respectively. 

5. Temporal variability 

The temporal variability of young low mass sources is often 
complex: the most common phenomena are magnetic flares 
with rapid rise and slower decays, superposed on an appar- 
ently constan t or, sometimes, ro tationally modulated emis- 
sion JFlaccomio et all 1200.4 IWolk et all l2005h . Other forms 
of variability, with less clear physical origin, are however often 
observed. 

We first investigated X-ray variability using the non-b inned 
one -s ample Kolmogorov-Smirnov (KS) test JPress et all 
Il992l) . This test compares the distribution of photon arrival 
times with that expected for a constant source. The test was 
applied to photons in the source extraction regions, which also 

16 The use of the J-band is justified because (i) in the presence of 
disk excesses the J-band is the most representative of the photospheric 
emission and (ii) the mass ranges in which is the mass-luminosity re- 
lationship is degenerate are narrower than for similar relationship with 
the H and K bands. 



contain background photons. Given that the background was 
found to be low and constant (§2.1), the results, i.e. the con- 
fidence with which we can reject the hypothesis that the flux 
was constant during our observation, can be attributed to the 
source photons. Tabled column 17, reports the logarithm of 
the KS-test significance with values < -4 truncated at that value: 
sources with log(Pics) <-3.0 can be considered almost defini- 
tively variable as we expect at most one of the 1003 sources 
(i.e. 0.1%) to be erroneously classified as variable. Eighty-five 
X-ray sources («8.5% of the total) fall in this category. Sources 
with -2.0 < log(PKs) < -3.0 can be considered as likely variable, 
but up to 9 such sources (on average) might actually be con- 
stant. Forty-nine Cygnus OB2 sources fall in this category. 

These numbers of sources in which variability is detected 
are however a lower limit to the total number of variable 
sources in the region for several reasons: i) most of the ob- 
served variability is in the form of flares, i.e. events that are 
shorter than our observation and with duty-cycle that ma y be 
on the other hand considerably longer IWolk et all 120051) : ii) 
the sensitivity of statistical tests to time variability of a given 
relative amplitude depends critically on photon statistic. This 
is illustrated in Figure|5] where we plot the fraction of variable 
sources as a function of source counts: the clear correlation be- 
tween the two quantities is most likely due to this statistical 
bias even though we cannot exclude a real dependence. 
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Fig. 5. Fraction of variable sources as a function of source 
counts. Horizontal bars indicate the range of source counts. 
Filled circles indicate, for each counts bin, the fraction of 
sources with log(Pics)<-3.0 (85 sources in total), while the up- 
per tip of the vertical bars indicate the fraction with log(Pics)<- 
2.0 (134 sources in total). The continuous and dashed horizon- 
tal lines indicate the overall (average) variability fractions ob- 
tained for the two confidence levels. 

Next, we extracted binned light-curves for each of our 
Cygnus OB2 X-ray sources adopting a bin length of 600 sec- 
onds, a compromise between bins that are long enough to reach 
a good signal-to-noise ratio per bin for most sources and suffi- 
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Fig. 6. Light curves (in the 0.5-8.0 keV band) for 20 sources with variable emission during our 97.7 ksec Chandra observation. 
The bin length is 600 seconds. Source numbers are given in the upper-left corner of each panel, followed by an 'f for sources 
classified as flaring. A detailed study of X-ray properties in flaring sources and the flare frequency in the Cygnus OB2 region will 
be presented in a forthcoming paper (Albacete-Colombo et al. in prep.) 
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Fig. 7. Light curves in the 0.5-8.0 keV band for the two variable early -type stars: Left: Source Id. #60, known as a B5Ie, has been 
recently classified a B5+0.5 Ia+ star by Klochkova & Chentsov ( 2004J). X-ray variabilit y of this star was unknown til l now. The 
time bin-size is 900 sec. Right: Source Id. #979 is probably a B1.5 V+'?' binary system (M assev & Thompson! 1 1 99 ll) . The time 
bin-size here is 1300 sec. Source number and results of the KS test are given in the legend. 



ciently short to resolve the decay phase of typical flares. Since 
the background of our observation is both small (negligible for 
many sources) and constant, we did not apply any background 
subtraction. We inspected the 857 light-curves of source with 
>10 photons, finding flare-like events, qualitatively defined as 



a rapid rise and a slow decay, in 98 sources, i.e. as 9.8% of the 
cases. These sources are indicated by a f in the last column of 
Tabled Of these 98 sources, 66 (-65 %) have a log(P KS ) < - 
3.0, while 13 were classified as probably variable (-2.0 < log 
Pks^-3.0). The remaining 19 sources were not detected as 
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variable by the KS test (log P K s>-2.0) and may or may not be 
actually "flaring". Figure[6] shows light-curves for 20 variable 
sources, 14 of which classified as flaring. Some of these sources 
(e.g. #1, #676 and #911) experience "impulsive" flares with 
very quick rises and decay phases of only a a couple of hours. 
Others (e.g. #33, #172, and #260) show longer (2 to 10 hours) 
flares. In several instances a second impulsive event is visible 
during the exponential decay of a previous flare (e.g. sources 
#6, #52, #439, #600, #796, and #945). Other sources (e.g. #555, 
#834 and #863) have variable light-curves that bear little re- 
semblance to typical flares and are instead characterized by 
slow continuous rises or decays that might be explained by ro- 
tational modulati on of n on-homogeneously distributed plasma 
jFlaccomio et alll2005h . 

Finally, we note that, as expecte d if flares or i ginate from 
magnetic reconnection events dFavata & Micelal |2003), the 
median photon energies for flaring sources are generally higher 
than those of non-flaring sources: the distribution of median 
energies for the flaring source peaks at E x &2.6 keV, with a lcr 
dispersion of 0.4 keV, while for non-variable stars it peaks at 
E x &2.1 keV, with a lcr dispersion of 0.3 keV. A similar conclu- 
sion can be drawn from the spectral modeling presented in the 
next section ({j6j: flaring sources often require higher tempera- 
ture models than non-variable ones. 

The variability of massive O and early B-type stars is sig- 
nificantly different from that of low mass members. Among 
the 26 OB stars detected in our FOV, 24 are classified as non- 
variable (log(P K s)>-2.0), in spite of their higher than average 
statistics, having between ~40 and 15000 counts, with a me- 
dian of ~ 1 1 1 counts. A comparison with Fig.[5]shows that the 
variability fraction of OB stars (7.7%) is significantly lower 
than for the bulk of our sources with similar statistics. This 
finding agrees with the common view of X-ray emission from 
O stars, which is believed to be unrelated solar-like magnetic 
activity, and rather explained as the integrated emission from 
a large numb er of small shocks occurring in the strong wind s 
of these stars [Feldmeie maD.ll997tlOwocki & Coherlll999l) . 
Interestingly however, two of these sources, #60 and #979, 
are significantly variable, with log(P K s) values lower than -4. 
FigureQshows their light-curves. 

Source #60 (Cyg#12) is a B5 Ie star JHansonil20ol . The X- 
ray light-curve shows a roughly linear decay of the ACIS count 
rate from ^0.18 cnts" 1 t o »0.16 cnts 1 . The star was observed 
bv lWaldron et alJ £l998) who did not report variability during 
-125 ksec of non-continuous ROSAT PSPC observation. This 
kind of variability may be similar to the rotational m odulation 
observed on the 07 star 9 OriC (Gag ne et all 12005) and there 
attributed to the presence of a rigidly rotating magnetosphere 
(RRM), which may form in magnetic e arly-type stars with mis- 
aligned magnetic and rotation axis dOwocki et all |2005). We 
defer further discussion about the nature of the Cyg #12 to 37.11 

Source #979 is identified with star #646 in lHansonl J2003h 
and there classified as a B 1 .5 V + ? star, the question mark in- 
dicating the presence of an unresolved faint secondary. The 
observed variable X-ray emission is much fainter than in the 
previous case and might be due to the magnetic activity of the 
companion to the B1.5 primary, presumably a lower mass star. 



6. Spectral analysis 

In order to characterize the hot plasma responsible for the X- 
ray emission of Cygnus OB2 stars, and to estimate their intrin- 
sic X-ray luminosities, we analyzed the ACIS spectra of each of 
our 1003 sources. Reduced source and background spectra in 
the 0.5-8.0 keV band were produced with AE (see 33.21 1. along 
with individual "redistribution matrices files" (RMF) and "an- 
cillary response files" (ARF). For model fitting spectra were 
grouped so to have a specified number of events in each en- 
ergy bin. Grouping was tuned to the source statistics and we 
chose 2, 5, 7, 10, and 60 counts per channel for sources with 
net-counts in the following ranges: [0-40], [40-100], [100-200], 
[200-500], and [500-16000]. Spectral fitting of background- 
subt racted spec tra was finally performed with XSPEC vl2.0 
llArnaudi r2004) and our own shell and TCL scripts to automate 
the process as described in lFlaccomio et alJ J2006I) . Best fit pa- 
rameters for the chosen models were found by x 2 minimization. 

We fit our spectra assuming emission by a thermal 
plasma, in collisional ionization equilibrium, as modeled 
by the APEC code (Smith etal., 2001). Elemental abun- 
dances are not easily constrained with low-statistic spectra 
and were fixed at Z=0.3 Zp, with solar abundances taken 
from Anders & Grevessdlll989i) . The choice of sub-solar abun- 
dances is s uggested by several X-ray studies of s tar forming re- 
gions (e.g. lFeigelson et alll2002t|Preibiscbll2003h . Absorption 
was accounted for using the WAB S model, parametrized by 
the h ydrogen column density, Nh JMorrison & McCammorl 
1983). 

Except for X-ray sources associated with O- and early B- 
type stars, which are discussed separately in 37.11 we fit source 
spectra with one-temperature (IT) plasma models using an au- 
tomated procedure. In order to reduce the risk of finding a rela- 
tive minimum in the \ 2 spaces, our procedure chooses the best 
fit among several obtained starting from a grid of initial values 
of the model parameters: log(N H ) = 21.0, 21.7, 22.0, 22.4, 22.7 
and 23.0 cm 2 and kT = 0.5, 0.75, 1.0, 2.0, 5.0 keV. Best fit val- 
ues of log(A^n) < 20.8 crrT 2 were truncated at 20.8 because, in 
the 0.5-8.0 keV energy range, ACIS spectra are insensitive to 
lower column densities. For the same reason, 154 best fit val- 
ues of kT turned out to be >10 keV and were truncated at that 
value. 

As noted by many a uthors ("e.g. lGetman et alll2002i l2005b: 
iFlaccomio et all 120061) . the spectral fitting of low statistic 
ACIS sources is problematic because of a degeneracy between 
plasma temperature and absorption. The degeneracy also re- 
sults in a systematic bias: kT values are often underestimated 
by as much as ~50% while Nh values are overestimated. We 
investigated this issue with our data by considering the distri- 
butions of the best fit parameters for source in different count- 
statistic bins. Figure [8] shows the run of mean kT and Nh with 
source counts for spectra fitted with IT models. The system- 
atic decrement of Nh and the increment of kT with increasing 
source statistics are hardly explainable as physical effects and 
rather indicate that the spectral parameters obtained for sources 
with less than ~20 photons are ill-constrained. 

The uncertainty on the Nh is particularly serious as it 
implies large and systematic uncertainties on the absorption- 
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Table 3. Results of X-ray spectral fits: first rows. The complete 
version is available in the electronic edition of A&A. 
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Fig. 8. Median Nh and kT vs. source counts for sources fit with 
one-temperature models. Filled and empty squares indicate the 
average values of the Nh and kT in each considered count range 
with vertical scales given on the right and left-hand side, re- 
spectively. The extent of the count ranges is indicated by the 
horizontal error bars. Vertical error bars indicate the median of 
the lcr uncertainties on the two parameters. For this plot we 
only considered N H and kT values with formal relative errors 
smaller than 90%. The hatched area below 20 counts indicates 
the count-range in which spectral fits were discarded because 
of strong biases in the best fit parameters. 



corrected X-ray luminosities. In order to reduce the risk of er- 
roneous results, we discarded from the following analysis re- 
sults for source with <20 net photons, a total of 423 sources. 
Table|3presents the results of the automated IT spectral fits for 
554 sources (the results of the spectral fitting for the 26 sources 
associated with know OB stars will be presented in 37.11 Table 
@}. We list source numbers (column 1) from Tables 1 and 2, 
background subtracted counts in the spectra (2), the best fit hy- 
drogen column densities and their 1 <x errors (3), the plasma 
temperatures and their 1 <x errors (4), and the emission mea- 
sures (5). Column (6) and (7) give the reduced xl for the spec- 
tral fits and the relative degrees of freedom, respectively. 

The hydrogen column densities derived from X-ray spec- 
tral fitting depend of the interstellar material in the line of sight 
to the Cygnus OB 2 cloud and on the location of stars within 
the cloud material of Cygnus OB2 . Figure[5Ja) shows the dis- 
tribution of log Nh values for the 580 sources with more than 
20 counts. They appear to be normally distributed with a me- 
dian log N H -22.25 (cm 2 ) and a FWHM of ~ 0.21 dex. The 
shaded area in Figure|9ja) indicates an apparent excess of rel- 
atively unabsorbed sources (log Nh <21.8 cm' 2 ) with respect 
to the log-normal distribution of the bulk of the sources. These 
~23 sources are likely associated with foreground stars. In fact 
we noticed that most of these sources seem to be spatially dis- 
tributed uniformly in the FOV of the Cygnus OB2 region. 
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0.78 


15 


8 


33 


21.70±0.08 


9.99 


53.40 
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Another interesting observation is that the Nh distribution 
of the 106 X-ray sources without near-IR counterparts (dark 
gray histogram in Fig.|9} seems to be skewed toward higher 
values with respect to global distribution. For example, while 
only 11% of the X-ray sources with counterparts have Nh > 
22.5, this is true for 29 (27%) of the unidentified sources. 

According to the relationship between Nh and Ay, 
N H =2.2xl0 21 A v cirT 2 dRvteJ. fl996). the median N H of all 
detected sources with more than 20 counts corresponds to 
A v ~ 8.li 2 j 8 mag. If we only consider the sources with near-IR 
counterparts (and > than 20 counts) we obtain log Nh=22.23, 
which corresponds to A v =7.7 mag., i.e. in reasonable agree- 
ment with the value calculated from the NIR CMD (A v »7.0 
mag). 

Figure|5Jb) shows the distribution of plasma temperatures: 
it peaks at ~1.35 keV, has median ~2.4 keV and shows an 
extended hard tail. Variable sources (logPKs < -2) have 
harder spectra with median kT=3.3keV and flaring sources 
(see hatched histogram in Fig.^) are even harder with median 
kT=3.75 keV. 

Finally we derived unabsorbed X-ray luminosities for each 
of our sources. For O and B-type stars, discussed in detail be- 
low ( 37. 1> luminosities were derived from individual spectral 
fits. For the other sources, given the considerable uncertainties 
in the absorption estimates, especially at the faint end, we pre- 
ferred not to use the results of spectral fits individually. Rather, 
we computed a single count-rate to Lx conversion factor that 
takes into consideration the average intrinsic source spectrum 
and interstellar absorption. This conversion factor, computed as 
the median ratio between the individual unabsorbed Lx's (from 
the best fit spectral models) and the source count-rates is 7.92 
10 33 ergs. 
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Fig. 9. (a): Distributions of absorbing columns, Nh, with corresponding visual absorption scale at the top. The light gray area 
refers to sources with net-counts greater than 20. The maximum of the distribution is at log Nh ~ 22.25 crrT 2 with a Full Width 
Half Maximum (FWHM) of about 0.21 dex. The shaded area on the left indicates an excess of low-absorption sources with re- 
spect to the log-normal distribution that best describes the histogram (the continuous curve), possibly associated with foreground 
stars. The darker gray histogram refers to the subsample of sources without a 2MASS counterpart, (b): same as panel (a) for 
plasma temperatures (kT) with, in addition, the distribution for "flaring" sources (hatched histogram). The peak of the overall 
distribution is at ~ 1.35 keV. 



7. Results 

We now discuss the implications of our data for the understand- 
ing of X-ray emission from young stars and for the study of 
the Cygnus OB2 stellar population. We split the discussion ac- 
cording to stellar mass and, correspondingly, expected X-ray 
emission mechanism. We first discuss known O and early B 
(<B2) stars, i.e. massive stars (M > 12M Q ) expected to gen- 
erate X-rays in their powerful stellar winds. Next, using our 
mass estimates based on the J-band magnitude ( 84.2t . we de- 
fine and study the properties of two subsamples of intermediate 
and low-mass X-ray detected stars. We base the distinction be- 
tween the two groups on the model-predicted thickness of the 
convective layer, whose presence has been often found to be 
correlated with solar-like coronal activity. Stars with small or 
absent convective layers are instead believed to be unable to 
sustain coronal activity, although the matter is somewhat con- 
troversial. Adopting the 2Myr isochrone from the SDF models, 
we observe that the thickness of the convective layer relative to 
the stellar radius, R c , decreases rapidly, from ~5% to <1%, as 
the stellar mass increases from 3 to 3.3M . This rather sharp 
boundary however lies right in the 2.8-5.3M range in which 
our mass estimates are degenerate and therefore unconstrained 
(see 34.2> . In order to resolve this ambiguity, we define the 
intermediate-mass range as 2.8-10 M , resulting in the selec- 
tion of 57 stars 17 . Low mass stars are finally defined as stars 
with M<2.8M and are predicted to have substantial convective 

17 There are only 28 stars with masses in the 2.8-5.3M Q range and 
for which this interval formally corresponds to the uncertainty on the 
mass estimate. We therefore expect a small fraction of the 57 interme- 



envelopes (R t . >18%). 578 counterparts to X-ray sources fall in 
this category. Because of the sensitivity limit of 2MASS only 
25(4) sources of these stars have masses below O.5(O.4)M , 
which we may consider as our approximate detection limit. 
Our completeness limit is however probably higher, at ~1M , 
corresponding J=15.8, i.e. roughly the completeness limit of 
2MASS. An histogram of the logarithm of stellar masses also 
peaks at ~1M , confirming the above estimate given that the 
mass function usually increases down to lower masses. 

Figure ^3 shows the L^-mass scatter plot for 683 X-ray 
sources in all three mass ranges. We distinguish with different 
symbols stars of low- and intermediate- mass, as determined 
from the J-band magnitude, and massive stars with known O- 
and early B- spectral types 18 . Note that stars in the 2.8-5.3M 
and 0.2-0.4 M ranges are plotted with horizontal error bars 
spanning the entire ranges, because the J-band vs. mass relation 
was there found to be degenerate. In the plot we also indicate 
the position of 57 time-variable X-ray sources. They appear to 
be brighter than other sources of the same mass, probably be- 
cause variability is more readily detected in sources with high 
photon statistics (see |5}- 



diate mass stars to lie in the 2.8-3.0M Q range and to have a substantial 
model-predicted convective envelope. 

18 A number of X-ray detected stars with no spectral type appear 
to be of high mass but with lower X-ray luminosities with respect to 
those with spectral types. In the discussion of high mass stars, how- 
ever, we will only consider stars with spectral types. 
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Fig. 10. Log M/Mq vs. log L x for Cygnus OB2 sources with 
NIR counterparts (black circles and star-symbols) and for ONC 
sources observed by the COUP project (gray circles). Variable 
sources are indicated by empty squares, known O- and early 
B-type stars by star-symbols. 



7.1. High-mass stars 

The Cygnus OB2 region is one of the most ma ssive SFR in 
the G alaxy, with -2600+400 OB star members ( Knod lsedeit 
l2000t) . The presence of a variety of evolved stars sug- 

[ ests that the star formation process was non-coeval 
Massev & Thompson! Il99ll iHansonl I2003I) . In fact, out of 
a total of 20 O-type and 13 B-type stars lying in our 
Cygnus OB2 ACIS FOV, 7 and 3, respectively, are evolved (lu- 
minosity class: I-III). In this respect it is interesting to note that, 
while we detect X-ray emission from all O stars regardless of 
evolutionary status, we have only detected 6 B stars, among 
which all the evolved ones. The detection fraction of evolved 
B star is thus 100% while it is only 33% (3 our of 9 19 ) for un- 
evolved ones. 

The most widely accepted explanation for the X-ray emis- 
sion of single O and early B-type stars invokes multiple small- 
scale shoc ks in the inner layers of their radiation-driven stellar 
winds fe.g. lFeldmeier et allll997l) . Recent theoretical and ob- 
servational results however support, for several OB stars, an- 
other physical plasma heating model involving strong magnetic 



fields: the ma 


jnetically channeled wind shock (MCWS) model 


( Schulz et all 


2003; Owocki et al. 


2005; 


Gagne et al., 


2005). 



Moreover, binary systems in which both components are O- 
and early B-type stars, can produce intense thermal X-ray emis- 
sion from wind-wind interactions as well as non-thermal X-ray 
emission from Inverse Compton scattering. 

We modeled the spectra of O and B stars with one- or 
two-temperature thermal plasma emission models (apec), ab- 
sorbed by neutral interstellar and circumstellar material (wabs). 



One undetected B stars has unkown luminosity class. 



Because the dense stellar winds of evolved OB stars are af- 
fected by the strong stellar UV/EUV ionizing radiation, the ad- 
ditional absorption of X-rays by a partially ionized stellar wind 
should be also considered. For evolved stars we then decided 
to model the combined effect of ISM plus wind mate rial with a 
warm absorption model (wabs x absori, see e.g. Wal dron et all 
1998). For sources #60, #544, and #729, we found that the in- 
clusion of a warm absorber in the spectral model yields satis- 
factory fits of the soft part of the spectra (<1.2 keV), reduc- 
ing the xl from 2.1 to 1.4, from 2.7 to 1.3, and from 2.2 to 
1.2, for the three sources respectively. The spectra of other 
evolved sources did not instead support the presence of the 
warm absorber, although in many cases the spectra have too 
low statistics to be conclusive. Metal abundances were fixed at 
Z = 0.3 Z Q in fitting faint sources (< 100 ph.) and left as a free 
parameter for brighter sources. When the best fit abundance 
was lower than Z = 0. 1 Z Q we repeated the fit with Z fixed at 
this value. 

Figure^] shows three notable examples of spectral mod- 
eling and the best-fit parameter values for the 26 detected OB 
stars are presented in Table 0] We list: X-ray source numbers 
from Tabl e 1 (column 1) ? identification numbers and spectral 
types from Hanson (2003) (2 and 3), reduced^ 2 of the best fit 
(4), fraction of pile-up as reported by AE (as discussed 33.21 
if f(Pile-up)>2% we used an annular photon extraction region 
that excludes the peak of the PSF), the N# of the cold and 
hot absorption components from wabs and absori, respectively 
(6 and 7), plasma temperatures of the two components (8 and 
9), metal abundances (10). Column (11) finally gives the un- 
absorbed X-ray flux in the 0.5 - 8.0 keV energy band. 

The median absorption of the 26 OB stars is log A^h=22.19 
ctrT 2 , very similar but slightly lower than the median value of 
lower mass stars (22.23, see ^7.31 . It is tempting to attribute 
the lower absorption of O stars to the sweeping of the cloud 
material by strong winds. A two-sided KS test comparing the 
distribution of A^h for O and B stars with that of lower mass 
ones gives, however, a null result. 

Most of the O- and early B- type stars are well fit by a 
soft thermal emission model with kT ranging between 0.5 and 
0.7 keV (~5. 8-8.1 MK), in agreement with the predictions of 
the w ind shock model JLucv & WhiteL Il980t lOwocki et all 
Il988l) . The spectra of some sources, however, require a second 
(hard) thermal component, which is not predicted by the 
wind shock model. The presence of this component seems 
to depend on the stellar evolutionary status: out of the 10 
evolved O- and early B-type stars (spectral types B5 Ie, 04III, 
05If, 03If, 06 If +05. 5111(f)) 5 (50%) show evidence of a 
hard thermal emission with kT roughly between 1.3 and 3.0 
keV (14-35 MK), while among the 16 class V stars, the hot 
component is only required in 3 cases (i.e. 23%). Several 
physical mechanisms have been proposed to explain the hard 
emission observed in some of our OB stars: 



Magnetically Confined Wind Shocks (MCWS). X-ray 
emission from MCWS is predicted by simulations of 
radiation-driven winds in th e presence of a magnetic field 
Jud-Doula & OwockiL 120021) and it has most likely been 
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Notes: ( 1 ) Abundance values in parenthesis indicate cases for which abundances were fixed in the spectral fit (see text), (2) Absorption-corrected 
fluxes are calculated in the 0.5-8.0 keV energy range (3) Three bright sources, #60, #488, and #568, suffer of pile-up, and we extracted the 
spectra from annular rings (see 33, 2> that exclude the central parts of the PSFs. {") the 6.4+6.7 keV FeK„ complex is observed; (*) the Fei 6.4 
keV line is not detected, but the 6.7 keV seems to be present. 



obser ved in O-type stars (e.g. fl'OrionisC, iGagne et all 
120051) . The model predicts that the overall degree of plasma 
confinement is determined by a single dimensionless pa- 
rameter 77+ = -jf^ , where B eq is the equatorial dipole mag- 
netic intensity at the stellar surface, R* is the stellar radius, 
M is the mass-loss rate, and Voo is the terminal wind speed. 
Simulations with typical parameters for OB stars predict 
high post-shock temperatures of about 20 - 30 MK (1.7-2.6 
keV). Furthermore, in the case of oblique rotators, i.e. when 
the rotational and magnetic axes are misaligned, rotational 
modulation of the X-ray emission can be observed. This 
might be the case of our source B5 Ie (Id. #60 = Cyg#12) 
for which kT2=1.56 keV and the light-curve shows, during 
our observation, a gradual decay (c.f. Figs. 171 and ll2b . 
- Highly compressed wind shocks in evolved O-type stars or 
colliding wind s in th e case of early - type b inary systems 
JStevenset all Il992l iRaassen et all l2003h . Our source 
#568 (Cyg#8) is an 06If + 05.5111(f) supergiants-giant 
pair, whose X-ray spectrum indicates hot thermal plasma at 
kT2 = 2.64+0.07 keV. A colliding wind scenario is more- 
over suggested by the presence of the FeKc, complex 



(6.4+6.7 keV) in its spectrum 20 (c.f. Fig. ll2t . The other two 
O-type super-giants, sources #544 (03 If) and #625 (05 If, 
see Fig. I12> . although requiring a second thermal compo- 
nent, show rather soft spectra and no FeKa, (6.4+6.7 keV) 
complex. 

- Inverse Compton scattering of photospheric UV photons 
by relativistic particles that are Fermi-accelerated in shocks 
with in the chaotic s tellar winds of OB supergiants (c.f. 
IChen & White! fl99ll) . The actual production of the needed 
relativistic electrons by stellar wind shocks is, however, still 
in doubt. We can likely rule out this possibility for our 
Cygnus OB2 OB stars, given that all of them are well fit 
by thermal emission models in which the lines of several 
elements give a significant contribution. 



20 The temperatures reached by hydrodynamic shocks in the winds 
of single stars are usually not high enough to produce significant Fe K ff 
emission. In wide early-type binary systems, however, stellar winds 
collide with velocities close to and are thus heated to sufficiently 
high temperatures. In such cases, the cooler surrounding wind mate- 
rial, excited by the high energy radiation, can produce the 6.4 keV 
fluorescent Fe line. The observation of the Fe-compl ex can thus be a 
used as diagnostic for colliding-wind binaries iRaassen et alll2003T) . 
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One of the most intriguing and so far un-explained obser- 
vational results regarding the X-ray emission of O- and early 
B- type stars is the existence of a simple proportionality be- 
tween the X-ray luminosity (L x ), a nd the bolometric luminos- 
ity (Lboi) (Berghoefer e t all Il997l) . We explore this relation- 
ship for Cygnus OB2 O- and early B-type stars in Fig. 
L x values were calculated from our X-ray spectral analysis. 
Bolometric luminosities with relative uncertainties were esti- 
mated form the k nown spectral types and the calibrations given 
by Martins et ai1d2005b . 

The proportionality between L x and Lboi is retrieved, albeit 
with a large dispersion. However, when we distinguish between 
main sequence stars and evolved ones (i.e. giants and super- 
giants), we discover, for these latter, an apparently tight L x - 
Lboi relationship which is not a simple proportionality. In fact, 
for evolved stars, the xl with respect to the Lx=10~ 7 Lb o i re- 
lation is 2.61, while if we adopt the best-fit power law (with 
index F~1.8) the xl is reduced to ~1.3 . We test the statisti- 
cal significance of including this additional degree of freedom 
in the L x /Lhni rel ation. We used a F^-test as described in (e.g. 
lBevingtonlfl969T) . We found that F^-test ~ 0.6, which does not 
indicate a significant improvement of the quality of the fit (F^ 
should be greater than 1.1), even when adopting a significance 
level as large as 0. 1 . However, this result is affected by the small 
number of evolved O and B stars in our sample. A similar anal- 
ysis with more stars from other regions will be needed in the 
future to confirm, or disprove, this tentative result. 

A more detailed study of the physical origins of the X-ray 
emission from OB stars and on its dependence on stellar wind 
parameters will be presented in forthcoming paper. 



Linear fit for evolved stats: — 6_ 
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log L bo , 

Fig. 11. X-ray versus bolometric luminosity for main-sequence 
and evolved OB-type stars (open and filled symbols, respec- 
tively) in the Cygnus OB2 region. The dotted and continuous 
lines indicate log L x /Lb i= -6, -7, and -8. Evolved stars (i.e. 
giants and super-giants) seems to follow a power-law relation 
with index ~1.8 (dashed line). 



7.1 .1 . Specific sources 

We now discuss in some detail the results of our analysis for 
three noteworthy OB stars in the ACIS FOV that have been the 
subject of previous investigations. 

- Cygnus OB2 #8 (ACIS source # 568). It has recently been 
classified as an 06If + 05.5111(f) binary with a 21. 9d pe- 
riod a strong evidence of phase-locked X-ray variability 
JPe Becker & Rauwl I2005I) . Our observation only covers 
~5% of the orbital period and we do not de tect variability. 
Cygnu s OB2 #8 has also been observed by IWaldron et alJ 
(2004) with the Chandra High Energy transmission Grating 
Spectrometer (HETGS). They fit the X-ray spectrum with 
a a 2T thermal model with Ti=3.98±1.67 MK (0.34+0.14 
keV) and T 2 13.88+1.80 MK (1.2+0.15 keV) and obtain 
an L x /Lboi ratio of ~1.810~ 7 . The spectral parameters are 
in rough agreement with those obtained here (see FigurelT2l 
and results presented in Table|4}, but we find a somewhat 
higher L x /Lb i, ~4.0xl0~ 7 . The discrepancy is however 
due exclu sively to the d ifferent adopted bolometric lumi- 
nosities: IWaldron et ail J2004 use log L bo i = 39.79 from 
Biegin g et al.l (Il989h- while we her e adopt log(L bo i) = 39.4 
obtained from lHerrero et al] ll2002ll . 

- Cygnus OB2 #9 (ACIS source #488). An 05If super gi- 
ant, it is the strongest and most variable n on-thermal ra- 
dio emitter in Cygnus OB2 dVan LooH2005h . In X-rays it is 
however fainter than other evolved stars. Cygnus OB2 #9 
has not been indicated as a binary in the literature. 
However, periodic radio variability has been recently ob- 
served with a period of ~ 2 . 35 yr (Blomme, private com- 
munication). Wald ron et al.l (fl998) did not detect variabil- 
ity in EINSTEIN and ROSATX-my data. lRauw etaD (f2005) 
successfully fit the XMM-Newton EPIC X-ray spectrum 
with an absorbed 2-T thermal model (kT; = 0.63+0.03 and 
kT 2 ~ 2.4 keV). These results are roughly consistent with 
ours (c.f. Table@}. 

- Cygnus OB2 #12 (ACIS source # 60), has been variously 
classified and its spectral type is probably variable. Most 
investigations howe ver agree with an extremely luminous 
B supergiants B5 Ie (Hanson, 2003), possibly related to the 
Luminous Blue Variable (LBV) phenomenon. Recent two- 
dimensional spectral classification suggests a B5+/-0.5Ia + 
spectral type (Klo chkova & Chentsovl 120041) . In our near- 
IR color-color diagram (Fig. [3}, Cygnus OB2 #12 appears 
to be highly absorbed and to be a ffected by a near-IR ex - 
cess typical of Herb ig AeBe stars llLada& Adamsi lT992). 
iKlochkova & Chentsovl J2004 present evidence of a line 
radial velocity gradient that can be interpreted as an indi- 
cation of matter infall. These findings suggest the existence 
of an accretion disk around the star. 

The XMM-Ne wton EPIC spectrum of Cygnus OB2 #12 was 
also fitted by iRauw et all d2005l) . who find temperatures 
kT, =0.73+0.16 and kT 2 = 1.8+0.4 keV, consistent with 
our own estimation (see Figure[Q]and Table0}. The tem- 
perature of the hard component is significantly higher than 
expected according to the wind shock emission model, 
given the low velocity of the stellar wind, ~ 150 kms" 1 . 
Emission from MCWS might better explain the spectrum 
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(see above). Such a scenario could also explain the ob- 
served X-ray variability (see |5} if the magnetosphere is 
tilted with r espect to the rotat ion axis as in the case of 
f3 1 Orionis C. lGaene et aTl l l2005l) . 



7.2. Intermediate-mass stars 

Intermediate-mass stars (2.8 < M/M < 10) are not expected 
to emit X-rays because, unlike O- and early B stars, they do not 
drive strong stellar winds nor possess outer convective zones 
that can sustain a dynamo mechanism such as the one that 
is ultimately held responsible for X-ray activity in low-mass 
stars. Several studies have however observed X-ray emission 
apparently associated with late B- and A-type stars. Although 
its origin has been often attributed to the coronal activity o f 
unresolved late-type companions (e.g. IStelzer et alll2006albl) . 
the matter is controversial. The Cygnus OB2 region has a large 
number of intermediate mass stars. It is therefore an ideal tar- 
get to check this hypothesis. We classify 66 ACIS sources as 
intermediate-mass. The Lx vs. mass plot in Fig. ^2 does not 
reveal any dramatic discontinuity at the boundaries with low 
and high mass stars. Intermediate mass stars, however, do not 
seem to follow the trend of increasing Lx with mass of lower 
mass stars, and seem to have Lx values that span the whole 
range spun by these latter. This observation is consistent with 
the companion scenario. 

As further supporting evidence for this conclusion we note 
that 8 of the X-ray sources associated with intermediate mass 
stars (about 12% of the sample) are variable, which is al- 
most the same fraction found for low-mass stars (see §6.3). 
Moreover, their X-ray spectra are well represented by isother- 
mal models with median kT=2.6keV, also in agreement with 
values found for low mass stars (see $1.31 

At the distance of Cygnus OB2 an offsets between IR and 
X-ray positions of 1", i.e. roughly the on-axis Chandra spatial 
resolution, corresponds to a projected separation of 1450 AU. 
X-ray source positions can however be determined with a preci- 
sion of ~0.1" (145 AU) or even better, depending on the source 
statistics and off-axis. If the intermediate-mass primary of a 
binary system is X-ray dark we should thus be able to distin- 
guish, at least statistically, whether the observed X-ray emis- 
sion originates from the low mass secondary in the case of wide 
binary systems, defined a s syste ms with separation > O.OOlpc 
(206 AU) (cf. lClose et all ll990). FigurefL3*l shows an example 
of an X-ray source associated with an intermediate-mass coun- 
terpart, but with a significant spatial offset. A comparison of 
the distributions of X-ray-NIR offsets for X-ray sources asso- 
ciated with low and intermediate mass counterparts is however 
inconclusive. We conclude that, if the companion hypothesis is 
correct, most X-ray emitting low-mass companions are not in 
wide systems. 

7.3. Low-mass stars 

The origin of X-ray activity in low-mass PMS stars has so 
far eluded full understanding. Although many aspects of MS 
coronal activity are also not understood, the X-ray activity 
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Fig. 13. Example of a possible misidentification of an X-ray 
source with a intermediate-mass near-IR star. We show the X- 
ray and K s -band images of ACIS source #216. The X-ray and 
2MASS positions differ by -2.5 arcsec (~1.8xl0~ 2 pc at the 
Cygnus OB2 distance). The X-ray source might be instead as- 
sociated with a faint near-IR source undetected in 2MASS. 

of MS stars correlates well with stellar rotational and con- 
vection properties, thus supporting the picture of a solar- 
like corona that is ultimately powered by a stellar dynamo. 
For Pre-Main Sequence (PMS) stars, however, no such a 
correlation has been found and the picture appears some- 
what complicated by the presence of accretion disks, which 
seem to play an important role in the X-ray emission. On 
one hand plasma heated by shocks due to the accretion pro- 
cess is likely to be a positive source of so ft X-rays (cf . 
[ Kastner et all l2002l Istelzer & SchmittL 120041 IFavata et all 
2005 HFlaccomio et alll2006l) . On the other hand the presence 
of accretion disks seems to lower the overall activity level, al- 
though probably makin g the X- r ay em i ssion harder and more 
time variable (IFlaccomio etail l2003t IPreibisch et all l2005t 
IFlaccomio et alll2006l) ~ 

We have defined the low mass range as M<2.8M , based 
on the presence of a significant convection envelope. The X-ray 
spectra and luminosities of low mass stars in Cygnus OB2 are 
quite typical of T-Tauri stars in other regions: our spectral anal- 
ysis (|6j yields a median kT=2.38 keV and logNH=22.23. 
Figure^|indicates that, although the spread of L x at any given 
mass is of the order of 1 dex, a trend of increasing L x with 
mass is observed. The same figure shows for comparison the 
nearly complete sample of Orion Nebular Cluster (ONC) stars 
as observed by the Ch andra Orion Ultra-deep Project (COUP) 
jGetman et all l2005bl) . The sensitivity limits of the two ob- 
servations are quite different because of the larger distance 
of Cygnus OB2 (-1450 pc vs. 470 pc for the ONC), and the 
longer exposure time of the COUP observation (~838 vs. ~97 
ksec for our data). However, as already discussed, the most se- 
vere limit to the completeness of our sample of likely members 
with near-IR counterparts is likely to be the sensitivity limit 
of the 2MASS data. Note indeed that the Cygnus OB2 stars in 
Fig-E3 nave a well defined limit in mass, at ~ 0.4 - O.5M , 
rather than in X-ray luminosity. 

Figure [TUIshows that the luminosities of low mass stars in 
Cygnus OB2 and in the ONC are fully compatible. We note 
here that the numbers of stars in the l.O-2.8M mass range 
in the Cygnus OB 2 and ONC regions differ considerably. In 
Cygnus OB2 we count 334 such stars, about 10 times as many 
as the in a sample of 375 COUP-detected lightly absorbed 
(A v <2) ONC stars with mass estimates, in which we find 34 
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Fig. 12. ACIS-I spectra of three OB stars. The data were fitted with absorbed two temperatures optically thin thermal plasma 
models. Left: Cyg#12 (Src.Id. #60), showing a prominent FeK 6.7 keV emission line, indicative of hot thermal plasma. Center: 
Cyg#8 (Src.Id. #568), a colliding wind binary with a strong FeK a blend at 6.4+6.7 keV, arguing in favor of th e therma l nature 
of the hard part of the spectrum and consistent with the expected emission from a colliding wind region (Raa ssen et all [2003). 
Right: source #625 (05.5 If), a typical super-giant fit with a 1-T model with a soft temperature model: KT; =0.55+0.06 keV. 



stars in the same range JGetman et alll2005bl) . Now, assuming 
that the IMFs in the two regions are the same and that the con- 
sidered sample of lightly absorbed ONC stars is complete, we 
can estimate the total number of stars in our Cygnus OB2 ACIS 
FOV as 375/34x334-3700. Considering that our FOV contains 
only one fourth of the 106 OB stars known in Cygnus OB2, we 
can then scale this number by the number of OB stars to get 
the total expected Cygnus OB2 population: ~15000 stars. The 
total stellar mass of the cluster can be then estimated by mul- 
tiplying this number by the average stellar mass in the lightly 
absorbed COUP sample (0.84M o):M d ~ 1.2 x 10 4 M o , some- 
what lower than the estimate of Knodlseder (2000), who de- 
rived 4-lOxlO 4 M . 

Deeper and more extended X-ray and near-IR surveys of 
Cygnus OB2 will be invaluable for the study of young stars of 
all masses with un-precedented statistic. 

8. Summary and conclusions 

We have reported the results of a deep Chandra X-ray observa- 
tion pointed toward the ~2 Myr old star forming region Cygnus 
OB2. Source detection was performed using PWDdetect, a 
wavelet-based algorithm, supplemented by visual inspection, 
identifying 1003 X-ray sources in the 17' x 17' ACIS-I field of 
view. Data extraction was performed using the semi-automated 
IDL-based ACIS Extract package, which is well suited to the 
analysis of observations of crowded fields such as ours. 

The X-ray source list was cross-identified with optical and 
near-IR (2MASS) catalogs: 26 X-ray sources were identified 
with optically characterized OB members of Cygnus OB2 (out 
of the 33 lying in the FOV) and 775 with 2MASS sources. 
Among these latter sources almost all are believed to be 
Cygnus OB2 members. About ~80 X-ray sources without opti- 
cal/NIR counterpart are estimated to be of extragalactic nature 
(AGNs) while the remaining X-ray sources with no counter- 
part are likely associated with members that are fainter than 
the 2MASS completeness limit. 

In order to characterize the previously unidentified likely 
cluster members with NIR counterpart, we placed them 
on NIR color-magnitude (K s vs. H-K s ) and color-color (H- 
K s vs. J-H) diagrams. A first estimate of interstellar extinction 
was obtained adopting a 2 Myr isochrone for the low- and 



intermediate-mass stars and assuming that O- and early B-type 
stars lie on the MS. We find a median visual absorption for OB 
stars of A v ~5.6 mag, while low mass likely members seem 
to be slightly more absorbed, A v ~7.0 mag. We also use the 
2 Myr isochrone and the J magnitude to estimate masses of 
likely members assuming that they share the same distance 
and absorption. Our sample of X-ray selected members with 
near-IR counterparts reaches down to M=O.4-O.5M , and is 
likely complete down to ~1M Q . From the H-K s vs.J-H dia- 
gram we estimate the fraction of low-mass stars with NIR ex- 
cesses (23/519), finding it to be quite low with respect to the 
1 Myr old ONC region indicating that either (i) the disk frac- 
tion decreases steeply with stellar age, or (ii) the intense photo- 
evaporating UV-field due to the large population of O stars in 
the region has reduced the disk lifetime. 

At least 85 or 134 X-ray sources, i.e. -8.5% and 13% of the 
total, were found to be variable within our observation with a 
confidence level of 99.9% or 99%, respectively. The fraction of 
variable sources increases with increasing source counts, likely 
as a results of a statistic-induced bias. The light-curves of 97 
sources indicate an impulsive, flare-like, behavior while the rest 
show more gradual variations of the X-ray emission. Only two 
of the 24 detected O- and early B-type stars are detected as 
variable during our 97 ksec observation, in spite of the high 
statistics of their light-curves. These exceptions are the B5 Ie 
star Cygnus#12 (our source #60), showing a rather linear decay 
of the count rate during the observation, and a B1.5 V+? star 
(#979), which also shows a slow variability, maybe related to 
rotational modulation. In this latter case the detected emission 
may originate from the primary B 1 .5 component and/or from 
its presumably low er mass compa nion indicated by the spectral 
type published by Hansor3 (l2003h . 

We modeled the ACIS X-ray spectra of sources with more 
than 20 photons assuming absorbed thermal emission models. 
Median log(Ne) and kT values of the sources are 22.25 (cm" 2 ) 
and 2.36 keV, respectively. Variable sources, Pks < 0.01 have 
harder spectra (median kT: 3.3 keV) respect to non-variable 
ones (median kT:2.1 keV) and flaring sources are even harder 
(median kT:3.8). Sources associated with O- and early B-type 
stars are instead quite soft (median kT: 0.75 keV). Absorption 
corrected X-ray luminosities of OB stars were calculated from 
the best fit spectral models. For the other, typically fainter, 
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sources, we adopted a single count-rate to flux conversion fac- 
tor computed from the results of individual spectral fits. Low 
mass stars have L x 's ranging between 10 30 and 10 31 ergs -1 
(median 2.5 x 10 30 ). Variable low mass stars are on average 
0.4 dex brighter, probably because of a selection effect. These 
X-ray luminosities are consistent with those of similar mass 
stars in the slightly younger (lMyr) ONC region. O and B- 
type stars are the most luminous, with L x = 2.5xl0 30 -6.3xl0 33 
ergs -1 . Their X-ray and bolometric luminosities are in rough 
agreement with the relation L x /Lboi = 10~ 7 , albeit with an or- 
der of magnitude dispersion. The X-ray luminosities of the 10 
evolved OB stars in our sample however, seem to obey a bet- 
ter defined and steeper power-law relation with Lb i: L x ocL bo 8 j). 
Further studies of a larger number of evolved OB -type stars are 
needed to confirm of disprove this dependence. 
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